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A fluid interpretation of Cardassian expansion is developed. Here, the 
Friedmann equation takes the form = g{pM) where p]\j contains only 
matter and radiation (no vacuum). The function g{pM) returns to the 
usual STTpM/i^nipi) during the early history of the universe, but takes a 
different form that drives an accelerated expansion after a redshift z ~ 1. 
One possible interpretation of this function (and of the right hand side of 
Einstein's equations) is that it describes a fluid with total energy density 
Ptot = ~8^9{Pm) = Pm + Pk containing not only matter density (mass times 
number density) but also interaction terms px- These interaction terms give 
rise to an effective negative pressure which drives cosmological acceleration. 
These interactions may be due to interacting dark matter, e.g. with a fifth 
force between particles F ~ r"~^. Such interactions may be intrinsically four 
dimensional or may result from higher dimensional physics. A fully relativis- 
tic fluid model is developed here, with conservation of energy, momentum, 
and particle number. A modified Poisson's equation is derived. A study of 
fluctuations in the early universe is presented, although a fully relativistic 
treatment of the perturbations including gauge choice is as yet incomplete. 
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I. INTRODUCTION 



Recent observations of Type lA Supernovae [1,2], as well as concordance with other ob- 
servations, including the microwave background [3] and galaxy power spectra [4], indicate 
that the universe is flat and accelerating. Many authors have explored possible explana- 
tions for the acceleration: a cosmological constant, time-dependent vacuum energy such as 
quintessence [5-10], and gravitational leakage into extra dimensions [11]. 

Recently, Freese and Lewis [12] (Paper I) proposed an explanation for the accelera- 
tion which involves only matter and radiation, invoking no vacuum energy or cosmological 
constant whatsoever. In their model, called Cardassian, the universe has a flat geome- 
try as required by measurements of the cosmic background radiation [3] and yet consists 
only of matter and radiation. The Priedmann equation is modified from its usual form, 

= -^PM^ to 
pi 

H^-giPM), (1) 

where H = a/ a is the Hubble constant (as a function of time), a is the scale factor of 
the universe, and the energy density pm contains only ordinary matter and radiation. The 
function g^pu) reduces to -^rpu in the early universe, so that Eq. (1) reduces to the ordi- 
nary Friedmann equation during early epochs including primordial nucleosynthesis. Only at 
redshifts z < 0{1) does the function g{pM) differ from the ordinary Friedmann Robertson 
Walker (FRW) case; during these late epochs, g{pM) gives rise to accelerated expansion. In 

Paper I, the specific form of g{p) that was considered was Power Law Cardassian, 

Stt 



= TT^PM + Bpl,, (2) 



with 

n < 2/3. (3) 

The second term only becomes important once z < 0{1), at which point it dominates 
the equation and causes the universe to accelerate. Other possible functions g{p) [13] are 
discussed further below. 



There remains the question of the fundamental origin of these modifications to the Pried- 
mann equation. There is no unique four-dimensional or even higher-dimensional theory that 
gives Cardassian evolution. We consider two different motivations for these modifications: 

1) These functions may arise from fundamental theories of gravity in higher dimensions, 
as was discussed in [12]. Chung and Preese [14] showed that, generically, the Priedmann 
equations are modified as a consequence of embedding our universe as a three-dimensional 
surface (3-brane) in higher dimensions. 

2) Alternatively these functions may arise in a purely four-dimensional theory in which the 
modified right hand side of the Friedmann equation is due to an extra contribution to the 
total energy density. The right hand side is treated as a single fluid, with an extra contribu- 
tion to the energy-momentum tensor in (ordinary four dimensional) Einstein's equations. 
The two motivations may or may not be linked, in that the fluid interpretation may be 
intrinsically four-dimensional, or it may be an effective description of higher dimensional 
physics. 

In this paper, we restrict our discussion to four dimensions, and treat the right hand side 
of Einstein's equations as a single fluid. We consider models with an extra energy density 
associated with matter that contributes in such a way as to drive acceleration. This extra 
energy density may be intrinsically four dimensional or may serve as an effective description 
of higher dimensional physics. We take the total energy density of the matter 

Ptot = -^9{p) ^Pm + Pk (4) 

(plus possible internal thermal energy which is unimportant on cosmological scales) to con- 
tain not only the ordinary mass density pu (mass times number density) but also an addi- 
tional contribution px- For example, in Eq. (2), 

PK = (5) 

Given this total energy density, we can now compute the accompanying pressure, and find 
that the Cardassian contribution has a negative pressure, px < 0. This negative pressure is 



responsible for the universe's acceleration. In fact, one can obtain any negative equation of 
state Wk — Pk/pk < 0, including wk < — 1. 

The fluid approach has several advantages: (i) it is fully relativistic, (ii) it allows for 
the conservation of energy and momentum as well as of particle number, (iii) it admits a 
sensible weak-field limit which leads to a modified Poisson's equation, and (iv) it permits 
the study of fluctuations in the early universe, the study of effects on the cosmic microwave 
background anisotropics, and other observables. 

The primary purpose of this paper is to examine this fluid approach. However, we briefly 
speculate on a possible origin for this extra term px in the energy density. It may arise from 
(dark) matter self-interactions that contribute a negative pressure, for example through a 
long-range confining force which may be of gravitational origin or may be a fifth force. This 
self- interacting dark matter is different from any such component considered in the past, 
in that is has a negative rather than a positive pressure. We speculate on a form of the 
force between particles that may be responsible for such an interaction, F ~ r""^, although 
this Newtonian form must of course be modified on horizon scales. This description of a 
self-interacting dark fluid may be an effective description of a more fundamental theory. The 
fluid approach does not rely on the validity of such an interpretation of self-interacting dark 
matter, e.g., the interactions may be an effective description of higher dimensional physics. 

We begin by reviewing the idea of Cardassian expansion in Sect. II. We present a general 
fluid formulation in Sect. Ill, and then give specific examples in Sect. IV. In Sect. V we 
address the growth of density perturbations, and in Sect. VI we speculate on the possible 
origin of an interaction energy with negative pressure. 

II. REVIEW OF CARDASSIAN MODELS 

The general form of a Cardassian model was described in Eq. (1), in which a general 
function of matter density replaces the ordinary energy density in the Priedmann equation. 
The simplest version, the Power Law Cardassian model of Eq. (2), can equivalently be 
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written as 



SttG 



( 



J 



n-l 



Pm + 



Pm 



(6) 
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PCard 



For n < 2/3, the new term is negligible initially, and only comes to dominate at redshift 
z ~ 1 (once Pm ~ Pcard); once it dominates, it causes the universe to accelerate. We can 
consider the contribution of ordinary matter, with 



to this new term. Here, subscript refers to today. Once the new term dominates the right 
hand side of the equation, we have accelerated expansion. When the new term is so large 
that the ordinary first term can be neglected, the solution to Eq. (2) is 



so that the expansion is super luminal (accelerated) for n < 2/3. 

The Cardassian model also has the attractive feature that matter alone is sufficient to 
provide a flat geometry. The numerical value of the critical mass density for which the 
universe is flat can be modifled. For example, in Paper I it was shown that in the model of 
Eq. (2), the value of the critical mass density can be 0.3 of the usual value. Hence the matter 
density can have exactly this new critical value and satisfy all the observational constraints 
such as given by the baryon cluster fraction and the galaxy power spectrum. 

In a 'generalized Cardassian model,' other functions g{pM) of the matter (or radiation) 
density on the right hand side of the Friedmann equation can also drive an accelerated 
expansion in the recent past of the universe without affecting its early history [13]. Several 
of these alternative functions will be discussed below [see Eqs. (48) and (52) below]. 



Pm = PM,o{a/ao) 



(7) 
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(8) 



a oc 



III. BASIC EQUATIONS 



A. Perfect Fluid 



We use the ordinary four-dimensional Einstein's equations 
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G^, = 8nGT^,. (9) 

On the right hand side, we take as our ansatz that the energy-momentum tensor is made 
only of matter and radiation, and has the perfect fluid form, 

T^"" ^pg^'' + {p + p)uV, (10) 

where p is the pressure, p is the total energy density of the matter and radiation, and u'^ 
is the fluid four-velocity. Here the total energy density for matter includes not only the 
mass density pm (mass times number density) but also any interactions or additional terms 
(corresponding to the new terms on the right hand side of Eq. (1)). In general, p and p 
are functions of the mass density pu and of thermodynamic variables. We assume that at 
recent times matter is non-relativistic, i.e. the typical speeds involved are much smaller than 
the speed of light, but we do not assume that p <^ p. 

As our ansatz, we take the total energy density of the fluid during the matter dominated 
era to arise from the sum of three terms: 

P = Pm + Pinternal + PK: (H) 

P = Pm+Pk, (12) 

where 

Pm = muM (13) 

is the ordinary matter density of some particle of mass ui and number density n^, Pintcrnai 
and Pm are the ordinary internal energy density and pressure of matter (for example, for 
an ideal monoatomic gas at temperature Tm, Pinternai — ^^mTm and pm — tt^mTm), and px 
and pk are extra (Cardassian) contributions to energy and pressure. With regard to the 
gross properties of the universe (such as its expansion) , we can ignore Pmternai compared to 
Pm and pm compared to px for nonrelativistic matter; however, in the context of galaxies, 
Pinternai a^d Pm cau bc important. 
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Eq. (11) is quite general. Even in those cases where the total energy density is not a 
simple linear sum of terms (see Sections IIB and IIC below), one can always write Eq. (11) 
in this fashion. 

In Cardassian models we assume there is no vacuum contribution to the energy density. 
We also assume that the new (Cardassian) contributions pK and pk are only functions of 
Pm] i-e., the Cardassian fluid is barotropic. For example, in Paper I, we took 

PK = hpl,, (14) 

where 

— pLb (15) 
Stt 

with B as in Eq. (2). 

Note that, throughout the rest of this paper, we will focus on the matter dominated era 
and hence concentrate on the matter contribution to the fluid (rather than the radiation). 
With the ansatz of Eq. (11) and (12), Eq. (1) can be rewritten 

H"^ = SttGp/S. (16) 

Hence we recover the ordinary FRW equations, but with a modified energy density on the 
right hand side. 

One can think of this total energy density as including the effects of an interaction term: 
perhaps this term is simply an effective term on large scales (possibly arising from extra 
dimensions) describing the expansion of the universe; or perhaps this term is due to the 
interaction energy of the dark matter. A possible origin of an interaction energy with a 
negative pressure is described in Sect. VI. We iterate again that the fiuid approach is only 
one of the many ways that Cardassian expansion of Eq. (1) could result. 



B. Conservation laws 



The Bianchi identities guarantee the conservation of energy and momentum. 
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T"".^ = 0. (17) 

One can follow the evolution of the energy density along each fluid world line using the gen- 
eral relativistic fluid flow equations. In a comoving frame, energy-momentum conservation 
gives the (fully relativistic) energy conservation and Euler equations [15,16] 

p=-<^(p + p), (18) 

where the dot denotes a derivative with respect to comoving time and the tensor h^^ — 
9ii.v — UfjUi, projects onto comoving hypersurfaces. 

We impose in addition that mass (or equivalently particle number) is conserved, 

{PMU%^ = 0. (20) 

This will give us the usual dependence of the matter mass density on the scale factor of the 
universe. 

C. Thermodynamics 

General thermodynamic relations connect p and p. In particular, consider the first law of 
thermodynamics, TdisV^ — d{pV^ +pdV. In an adiabatic expansion, d{sV) — 0. Together 
with mass conservation, d{pMV) = 0, this equation leads to 

"^'"{i^h'' '''' 

which allows p to be computed from an expression for p. It also allows the speed of sound 
to be written as 

-^Sl-lSl/lfe). 

We can define w — p/p, which in general will not be constant. 

Using the ansatz given in Eqs. (11) and (12) in Eq. (21), we can now relate the Cardassian 
contributions to energy and pressure via 



Pk = PM ( 1^ J - PK- (23) 

Again we can define wk = Pk/ Pki which in general will differ from w = p/ p and will not be 
a constant. 

D. Newtonian limit 

Now we obtain the basic equations in the Newtonian limit. In Minkowski space, we 
write u°' = 7(1, v) and the metric rj""^ = diag(— 1, 1, 1, 1), where 7 = — and v is the 

3-dimensional fluid velocity. Then, T^,^ = diag(— p,p,p,p). 

We can obtain the gravitational field equations. In a weak static field produced by a 
nonrelativistic mass density, the time-time component of the metric tensor is approximately 
given by Qoq ~ — (1 + 20). Here is the Newtonian potential for the gravitational field 
g = — V0. The Ricci scalar R ~ Gqo — 2V^0. Poisson's equation follows from R — SnGT^^/j, 
as 

V^(l) = 47rG{p + 3p). (24) 

We also have 

Vx^=0. (25) 

In the same weak field limit, and for a fluid moving non-relativistically, the energy 
conservation and Euler's equations can be found following Ref. [17], but without assuming 
p < p. Prom T'^^.p = we find 

^ + V-[{p + p)v\^0. (26) 
Prom T'^^.p — 0, we find Euler's equation, 

|- + (ff.V)ff=-4^-W. (27) 

ot P + P 

Notice that we do not assume p <^ p in the right hand side of Eqs. (24-27). 



With the additional constraint of particle number conservation, we have the ordinary 
continuity equation for matter, 

^ + V-M = 0. (28) 

Using this equation, and in those situations where the internal energy and pressure can 

be neglected (e.g. when considering the overall expansion of the universe), we can rewrite 
Eq. (26) in the following way: 

v-^pk+ \Pk+Pk - PM^]^ -^=0. (29) 

V OpM ) 

Note that Eq. (29) in the homogeneous background of our universe reproduces Eq. (23), 
here in the Newtonian limit. 

Hence our basic nonrelativistic equations are Eqs. (24-28). 

E. Friedmann-Robertson- Walker models 

As mentioned above, Friedmann-Robertson- Walker cosmological models take the usual 
form [17]. The scale factor of the universe a obeys the equation 

a 4:71 G , „ , 

- = ^P + 3p 30 

a 3 



and the Priedmann equation 



8nG k 



where k — 0,±1 fixes the curvature of the spatial sections and H — a/ a is the Hubble pa- 
rameter. As mentioned previously, we consider only a flat universe with A; = 0, as motivated 
by microwave background data [3]. Of course, the total energy density here contains new 
terms and is given by Eq. (11). 

Cardassian expansion requires that the universe be accelerating today. Prom Eq. (30), 
at matter densities of the order of the matter density in the universe today, pm ~ Pmq, we 
require po = p(pmo) and po = p{pmo) to satisfy 
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po + 3po < (32) 

so as to have an accelerating universe. With given in Eq. (23), one can see that ac- 
celeration results if pu — '^Pk + ^Pm§^ < 0. In the limit where px ^ Pm, one can see 
that acceleration results as long as pk goes to zero faster than (for px = bp']^^, this 
requirement becomes n < 2/3 as stated previously). 

The energy conservation equation in an FRW model is 

p^-3H{p + p). (33) 

Particle number conservation gives 

Pm = -SHpM, (34) 
from which we obtain the usual scaling pM ~ ci~^ that was used in Eq. (7). 



IV. EXAMPLES 

In this section we discuss three different examples of barotropic Cardassian models, where 
the new contribution to the energy density is a function only of the mass density (mass 
times number density), px = Pk{pm)- In another paper we examine the consequences for 
supernova data of these three models [18]. 



A. Original Power Law Cardassian Model 

In the original Cardassian model of Ref. [12], 

PK = bpl, = PM (^) " , (35) 
\ Pm J 

with n < 2/3 as in Eqs. (2) and (6). The pressure associated with this model in the fluid 
approach follows from Eq. (23) as 

Pm J 
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Notice that 



Pk = -(1 - n)pK- 



(37) 



This model therefore has a constant negative wk — Pk/ Pk — — (1 — n). Then 



, / PCard\ 

Pk + Pk pm- 

\ PM 



(38) 



The speed of sound in this model 



n{l — n) 



n 



V PCard / 



1-n 



(39) 



is not guaranteed to be positive. So this model should be considered as an effective descrip- 
tion at scales where > 0. 

1. Basic Equations in the Newtonian Limit for Power Law Cardassian Model 
For the Cardassian model, the basic Newtonian equations become: 



dv -t Vp ^ 

— + [vV)v^ — + g, 

at P + P 



(40) 



V • VpK + 



\ PM ) 



Pm 



(41) 



(42) 



V PM , 



1—n 



Pm 



(43) 



dp 



M 



dt 



+ V • {pMV} = 0. 



(44) 
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2. Problem on Galactic Scales 



The fluid approach to the model cannot be used on galactic scales. Prom Eq. (40) 
and (38), we see that 

^ VpAf + ^Pk ^ ^ .^gx 

dt pm [1 + n{pcaid/pMy~''] 

where we have dropped pm Pm and Pintemai ^ Pm in the denominator. In standard 
cosmology, one would have 

-77 = — — + 9, (46) 
at Pm 

where pm is quite small and gravity dominates. Here, the fluid version of the p" Cardassian 
cosmology has an additional term in both numerator and denominator. Since pcard 
Pm throughout most of the interiors of galaxies, we can ignore the second term in the 
denominator. However, the second term in the numerator can have drastic effects. Let us 
consider, e.g., the position of the Sun in the Milky Way. The local energy density of matter 
is roughly pM ~ lO^pcard- Hence we have pm ^ Pk, so that one might expect to be able to 
ignore any effects of Cardassian terms. However, these expectations are not true, because 
of the Cardassian pressure, which in this case dominates: \pk\ ^ Pm- Hence there is a new 
force that acts due to the second term in the numerator of Eq. (45), 

dv VpK . X /pCardV"" VpM ..„x 

— = --— = +n(l - n) —— ——. (47) 
new \ Pm J Pm 

This force destroys flat rotation curves (velocities tend to increase as one goes out to large 
radii in an unacceptable way). The argument in the preceding paragraph has been presented 
using nonrelativistic equations, and as such is incomplete. However, we have found that the 
result that a problematic new force arises on galactic scales remains true in a fully relativistic 
treatment based on the Oppenheimer-Volkov equation. 

The fluid p"' Cardassian case must therefore be thought of as an effective model, which 
applies only on cosmological scales. The examples in the following two subsections, on the 
contrary, can be treated as fluid models on all scales, including cosmological scales as well 
as galactic scales. 
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B. Polytropic Cardassian 



Another class of models has 



P — Pinternal + PCard 



(48) 



\PCard / 

with 7^ 0. This model can be used on all scales (see below), but it does not quite fit 
the criteria of Cardassian as defined in Eq. (1); at late times in the future of the universe, 
when pm <^ Pcard, this model becomes cosmological constant dominated with A = pcard- 
Phenomenologically, this energy density is very similar to a model that was derived earlier 
[11] motivated by gravitational leakage into extra dimensions. 
In this model, the pressure is 



P = PM- PCard 



1 + f^T 



1 



(49) 



\PCard / 

When the ordinary internal energy density Pinternai and the ordinary pressure pm can be 
neglected, this model obeys a polytropic equation of state 

P = -Pcard ( — ^ 1 , (50) 

\PCard / 

with negative pressure and negative polytropic index N — —1/q. For q > I, the speed of 
sound in this model is positive, 

= r^f^- (51) 



1 + 



Pm 

^PCard 

Thus this fluid model can be used on all scales. 



We must make sure that at the scales of galaxies and galaxy clusters the Cardassian 
pressure can be neglected compared to the ordinary pressure. At large matter densities, the 
Cardassian pressure is \pk\ — plf'^Pcard- ■'■n ^ galaxy or cluster with velocity dispersion a, 
the ordinary pressure is pm — Pm(^'^- We want \pk\Ipm — {pCaxd/ PmY /(^"^ ^ 1- Taking 
(7 ~ 300 km/s and assuming pk is unimportant out to ~ 100 kpc where pm ~ lO^pcard, 
this condition amounts to g ^ 3. It is remarkable that this value of q is compatible with 
supernova data [18]. 
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C. Modified polytropic Cardassian 



A Cardassian model that can be used on all scales is 



P — Pinternal + RM 



qv-i 



\ PM ) _ 



(52) 



For g = 1 this reduces to the original p" Cardassian model with n — 1 — v. The pressure 
follows as 

/ PCard \ 



P = Pm- I^PM 



PCard \ 
PM J 



(53) 



\ PM ) 

This model is interesting because the two parameters v and q are important on different 
scales. The parameter v sets the current value of w ~ — i/, and so can be chosen to fit the 
supernova data, while the parameter q governs the suppression of the Cardassian pressure 
at high densities, and can therefore be chosen not to interfere with galactic rotation curves 
and cluster dynamics. Concrete comparisons with data will be presented in [18]. 



V. PERTURBATIONS 
A. Newtonian theory 

Here we calculate the behavior of small fluctuations, using the nonrelativistic equations 
derived above in Eqs. (24 - 28). As discussed below, there are gauge choices that render 
the Newtonian theory inadequate, yet we can learn from it nonetheless. We will write down 
the general perturbation equations for any generalized Cardassian model, and then as an 
example will solve them in the form of the Cardassian equations (of Paper I) to illustrate 
the type of results that occur. 

For the zero-order solution (superscript '^^^) we take the simple spatially uniform solution 
with 

(54) 



(0) 

Pm 



Pm,o 



a{t) 



15 



a{t) 



7(0) 



47rG'(p(°) + 3pW) 



(55) 
(56) 



We now seek a perturbed solution (superscript ^^^) by adding to the zero-order solution 
the small perturbations , and g^^\ Following [17], we find that the hydrodynamic 

equations Eq. (40-36) then give, to first order in these perturbations, 

1 



^1) I ^ ^1) 



PM + PK 



Pm^'^ + 3-p'S+Pm^-v^''^0, 
a 



V7 .-r(i) _ 



V X 



0, 



d , „ , 
1 + ^ — [Pk + ^Pk) 

OpM 



(57) 
(58) 
(59) 
(60) 



We have dropped the superscript (0) from zero-order quantities. In deriving these equations 
we have neglected pintemai and pu with respect to pu- 

Wc take the perturbations to be adiabatic, so that the ordinary pressure perturbation is 
given by 

(61) 



(1) _ 2 (1) 

Pm — '^s Pm 



where Vg is the ordinary speed of sound. In the generalized Cardassian model with pressure 
and energy density related by Eq. (23), the perturbation equations require 

d'^pK 



(1) 



Pk — Pm-^-^ Pm 
^Pm 



(62) 



1. Newtonian perturbation equations for />" Cardassian case 



For the case of Paper I in Eq. (2), Eqs. (62) and (60) become 



V • = -AttG 1 + n{3n - 2)bp 



and 



(63) 



= bn{n - 1)pVV^^ 



1 Ji) 



(64) 
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B. Newtonian theory: non-expanding case 

First, let us consider the nonexpanding case with unperturbed velocity v — {a/a)r — 0. 
A combination of the perturbation equations (57-60) yields 

= ^Inew ^'^^M + 47rG'pM,new ^M- (65) 

This is the usual perturbation equation for 

Sm = ^, (66) 
Pm 



with Vg replaced by 



2 Pm / 2 , d pK \ 



and Pm replaced by 



d d 
PM,new = Pm — (P + 3p) ^ Pm + Pm ^ — {pk + ^Pk) ■ (68) 

OpM OPM 

The equations are spatially homogeneous, so we expect to find plane-wave solutions. We 

take 5m oc exp[i(/c • x — uit)] to find 

= k'^vl^^^ - 4:7lGpM,new (69) 

If u"^ < 0, then the perturbation is unstable to collapse. Hence we find the Jeans length 

Aj = — = VTT ' (70) 

Perturbations on scales larger than the Jeans length can collapse. 

In particular for the p" Cardassian case of Paper I in Eq. (2), we have 



l + nbpH^ ■ ^ ^ 

and 

PM,new = Pm [l - ^^(2 - 3n)p^^] . (72) 
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Then the Jeans length becomes 
Aj 



vr 



Vg — bn{l — n)p 



n-l 

M 



Gp 



M 



1 + nhp'^M^ 



n—l 



bn{2 - 2>n)p1t 



(73) 



The Jeans length for ordinary cold dark matter was tiny, and the modifications here do not 
change it substantially enough to make it interesting. 



C. Newtonian theory: expanding case 

Here we consider the expanding case with v — {d/a)f^ 0. We take the plane wave form 
for the solutions, 



PM{r,t) = PM{t)S{f,t), 



S{f, t) — S{t) exp 



ir ■ q 
a{t) 



(74) 
(75) 



After some algebra, we find 



•■ 2a ■ 

6 H 6 

a 



(76) 



where f^nc^ and Pm, new are given in Eqs. (67) and (68). 

For the Power Law p** Cardassian model of Eq. (2), which is an effective model on 
large scales, we can drop the (f term. We have solved the resulting equation numerically 
after changing independent variable from t to a{t). The solution for 5{a) (for modes with 
A ^ Aj) is given in Figure 1 for varying values of the Cardassian index n. In the figure, 
ocard IS defined to be the scale factor at which the Cardassian term starts to dominate in 
the Friedmann equation, i.e. when pM(acard) = PA'(acard)- Today, ao/acard is a factor of a 
few, but the figure extends to much higher values of a to show the behavior of the solutions. 
The vertical axis is normalized so that the solution becomes a/ocard at early times; namely, 
using the solution in Ref. [19] for the growth of matter fiuctuations in the radiation and 
matter dominated eras of a cosmologically fiat model, 

5{a) 2aeq 



(^init 3acard' 
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(77) 




10-2 0.1 1 10 102 IQ3 

^ A Card 

FIG. 1. Growth factor D{a) for matter density perturbations in the original px = hp%[ Cardas- 
sian model, as a function of the scale factor a. The normalization scale factor acard is defined by 
PM(ciCaxd) = PKicLCaxd)- Curves are labeled by the Cardassian index n. The curve with n = 1 
corresponds to the usual growth in a matter dominated universe, while n = corresponds to a 
cosmological constant. 
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where 5init is the initial value of the pertubation and Oeq is the scale factor at matter-radiation 
equality. One can see in Figure 1 that perturbation growth is suppressed once a > acard for 
< n < 1. As a reminder, n = 1 corresponds to ordinary matter whereas n = corresponds 
to a cosmological constant. To remind the reader, this figure corresponds to perturbation 
growth in an expanding Newtonian fluid in a box with the sides being pulled out. 

D. General- relativistic theory 

Since we do not neglect the pressure p with respect to the energy density p, the pertur- 
bations should in fact be studied in a general-relativistic theory. The problem arises of the 
choice of gauge. In order to define the energy density perturbation throughout spacetime, 
one needs to choose a set of hypersurfaces; i.e., one needs to choose a gauge. Two choices 
are common: 

(1) The comoving gauge corresponds to the set of comoving hypersurfaces, defined as those 
which are orthogonal to the comoving world lines, i.e., to the world lines which follow the 
flow of energy. 

(2) The synchronous gauge corresponds to the hypersurfaces which are orthogonal to 
geodesies^ . 

The value of density perturbations depends on the choice of gauge. However, typically, if 
one looks inside the horizon and well into the matter dominated era, the value of the density 
perturbation becomes the same in all gauges. 

In the case of Cardassian cosmology, we have the unusual circumstance that, even inside 
the horizon and well into the matter dominated era, the value of the density perturbation 
depends on the choice of gauge. 

The evolution of density perturbations in the universe using a fluid flow approach was 



^ There are in fact an infinity of synchronous gauges, but authors typically drop the so-called 
"gauge mode" solution so that the synchronous gauge becomes unique. 
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discussed in Ref. [16]. In the comoving gauge, with w — p/ p and the relation 5p/5p — 
(true for barotropic Cardassian models), we find that the fractional perturbation S — 5p/p 
of momentum k obeys the equation [16] 

H-^'d + [2 - 3{2w - cl)]H-^5 - -(1 - Qc] + 8w - 2>w^)5 = - f -^fI c^S. (78) 

2 yaH J 

The pertubation 5s in the synchronous gauge is related to S by [16] 

'■^ Sp 



d-ds = 3H(1 + w 



) / -^dt, (79) 



if one drops the "gauge mode." In ordinary cosmologies, where there is no Cardassian 
pressure term, the integral on the right hand side has a fixed value after matter domination, 
obtained by setting 5p — in the matter-dominated era. A few Hubble times into the 
matter-dominated era it becomes negligible compared to the time scales of interest so that 
the comoving and synchronous gauges become identical. However, the values of the density 
perturbations in the two gauges do not become equal in Cardassian cosmology, when the 
Cardassian pressure 6pK is present. This term contributes to the integral in Eq. (79) all 
the way to the present time, and the comoving and synchronous gauges are not identical 
even today. This creates a problem of interpretation for fiuctuations in the present universe, 
which must be addressed in future studies. 



VI. ORIGIN OF INTERACTION ENERGY WITH NEGATIVE PRESSURE 

Here we speculate on a possible origin for an interaction energy with a negative pressure. 
Dark matter particles may be subject to a new interparticle force which is long-range and 
confining, 

F(r) oc r"-^ (80) 

with a > 0. This force may be of gravitational origin or maybe a fifth force. 
To be more quantitative, let us write the new interparticle potential as 
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Uij = ^rg, (81) 

where rij is the distance between particles and A is a normahzation constant. The total new 
interaction energy of a system of N particles occupying a volume of radius R will be 

C/new ^ AN'^R^, (82) 

to within a numerical factor of order 1 dependent on the geometry. The total gravitational 
potential energy of the same system is, also within a factor of order unity, 

C/grav - (83) 

where M is the total mass of the system. To play a cosmological role at the present time, the 
new energy must be of the same order of the gravitational energy when R ~ Rh, the current 
size of the horizon. Imposing that C/new — t^grav at ~ Rh gives us the normalization 



A = -B^, (84) 

where m — M/N is the mass of a single particle. 

We can now find the magnitude of the new force on galactic scales. We have 

Thus the new force per unit mass on a particle of mass m at distance Rg from a system of 



N particles is of order 



|Vf/„.wl GM f R„Y-' 

" -52- ■ (86) 



new ^ ^ new 



m m Rjj \Rh 

Compared with the gravitational force 



~ ™ (87) 



this gives 



^ new / "-e 



-^grav ^ 
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^i^r'- (88) 



For Rg of the order of galactic scales, i.e. Rg <ti Rh the new force is negligible compared 
to the gravitational force. It is clear that for a > the new force is only important on 
very large scales. This Newtonian formulation must of course be modified at large distances 
because of the finite speed of hght and issues of causality. 

We want to comment on the equation of state of a system subject to long-range confining 
forces. If such a system reaches thermal equilibrium (that it does so in the presence of long- 
range confining forces is not at all clear), then simple statistical mechanics considerations 
based on the scaling of the partition function lead to the equation of state (see Appendix 
for details) 

P = - gP- (89) 

This is the equation of state of the force mediators. For example, the Coulomb force (al- 
though not confining) has a — —1 and equation of state p — p/3, which is that of photons. 
If q; = 1 (as in QCD) or a = 2, the mediators are strings and 2-dimensional objects, re- 
spectively, and their equations of state are p = — (l/3)p and p = — (2/3)p, which are those 
of a network of strings and of domain walls, respectively. Finally, one obtains the vacuum 
equation of state p — —p for a — 3. Notice that the Cardassian index n in the Power Law 
model is connected to the exponent a in the confining force law through a = 3(1 — n) 
(cfr. e.g. Eqs. (37) and (89)). That a confining force can give rise to an effective negative 
pressure is well-known in particle physics, where the MIT bag model is just such an effective 
description of quark confinement. The cosmological negative pressure may be an indication 
that our observable universe is in a big bag. This suggests that the dark energy may be the 
interaction energy associated to a long-range confining force. 

VII. CONCLUSIONS 

An interpretation of Cardassian expansion as an interacting dark matter fluid with neg- 
ative pressure is developed. The Cardassian term on the right hand side of the Friedmann 
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equation (and of Einstein's equations) is interpreted as an interaction term. So the total 
energy density contains not only the matter density (mass times number density) but also 
interaction terms. These interaction terms give rise to an effective negative pressure which 
drives cosmological acceleration. 

These interactions may be due to interacting dark matter, e.g. with a long-range con- 
fining force or a fifth force between particles. Alternatively, such interactions may be an 
effective description of higher dimensional physics. We have said that matter alone can 
be responsible for accelerated behavior. However, if the Cardassian behavior results from 
integrating out extra dimensions, then one may ask what behavior of the radii of the extra di- 
mensions is required. Similarly, if we follow a QCD bag or other description of self-interacting 
dark matter, one may wonder if an equivalent vacuum description can be constructed. Fur- 
ther work in search of a fundamental origin of Cardassian expansion must be studied to 
answer these questions in detail. 

A fully relativistic fluid model of Cardassian expansion has been developed, in which 
energy, momentum, and particle number are conserved, the modified Poisson's equations 
have been derived, and a preliminary study of density fluctuations in the early universe has 
been presented. 

One of the goals of this study is to allow predictions of various observables that will serve 
as tests of the model. The Cardassian model will have unique predictions, particularly due 
to the modifled Poisson's equations. For example, one can now calculate the effect on the 
Integrated Sachs Wolfe effect in the Cosmic Microwave Background [20] . In addition, one can 
now calculate the effect on cluster abundances as function of redshift. These predictions can 
then be tested against existing and upcoming measurements of these quantities. Comparison 
with existing and upcoming supernova data is being studied in another paper [18]. We 
reiterate that this fluid approach is only one of the ways that Cardassian expansion may 
result. 
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APPENDIX A: EQUATION OF STATE FOR A SYSTEM SUBJECT TO 

CONFINING FORCES 

We give here the details of the derivation of the equation of state p — — (a/3)p for the 
mediators of a confining interparticle potential U — Ar°'. Our derivation assumes that the 
gas is non-relativistic and in thermo dynamical equilibrium. We stress that it is not at all 
clear that a system of particles subject to long-range confining forces may reach thermal 
equilibrium. For our application to dark energy, it is not even clear that it could reach 
equilibrium on a time scale short compared to the age of the Universe. Therefore, the 
content of this Appendix may be of somewhat academic interest. Nevertheless, we present 
it for completeness. 

The partition function of a non-relativistc gas of particles subject to a confining 
interparticle potential U = Ar°' is 

where T is the temperature and V is the volume occupied by the system. 

If we rescale V — > and T — > A^T, and then change integration variables r — > Ar', 
p — > A"/^p', we can prove that the partition function scales as 

Z(AV, A"T) = A=^^+i"^Z(l^, T). (A2) 



p 



(Al) 
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Now the free (ideal gas) partition function 

^o(V^,r) = ^"(^)'" (A3) 



scales in the same way, 

Zo(AV, X'T) = A=^^+t'^^Zo(l/, T). (A4) 

It follows that 

Z{V,T) = Zo{V,T)Z,{TyV"), (A5) 
where Zi{x) is a function of the ratio T^/V"', which is invariant under the rescaling V — > X^V, 

Pressure, entropy, and energy density can then be computed from the free energy 

F = -T\nZ{V,T) = -T\nZo{V,T) - Tin Z^{T^ /V^) (A6) 

as 

— j =-NT + 2,NTf{T'/V'^). (A8) 

where 

The first terms on the right-hand sides of eqs. (A7) and (A8) correspond to the ideal gas. 
The second terms are the pressure p and the energy pV due to the confining forces. They 
obey the relation 

P=-fp- (AlO) 

Although we have not quantized the interaction, we draw on the analogy with electromag- 
netism described in the main text and call eq. (AlO) the equation of state for the force 
mediators. 
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